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ABSTRACT 

The task of teaching principles systems analysis'and its applications 
to students with no programming background can be greatly achieved using 
high level simulation language like CSMP/36O. Various stochastic situations 
represented by Monte-Carlo simulation can be designed in FORTRAN with some 
pre-programmed segments given to the students. Heavy load of computer 
bomework assignments contributed to a much better understanding and to the 
speed of learning and higher motivation. Example of such an ass-i gnment is 
provided. 
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COMPUTER ASSISTANCE IN TEACHING DYNAMIC-STOCHASTIC SYSTEMS ANALYSIS 

Hovav TalpazsV o 
Introduction 

The systems concept has come to play a critical role in contemporary 
science [2,4] and has become an indispensable component of not only engi- 
neering but also the behavioral studies, especially at the university 
graduate program level [3,6,7,8]. 

Though an old concept, systems analysis became a useful discipline 
only with the appearance of the electronic computer, which, in combination 
with systems methodology [2,7], made it possible to ^olve complex problems 
with wide areas of application. Yet, another decade was to pass before the 
computer became general ly aval lable as an educational tool for use in systems 
analysis cQurses. 

In Agj-icul tural Economics, the volume of computer use would probably 

rank reseajrch first, extension second and teaching third. It is beyond the 

1 

scope of this paper to investigate the reasons for this. A few of the reason; 
include physical i naccessab i 1 i ty of computer services, funding structures, 
lack of instructors' experience, and conservatism in teaching methods as 
demonstrated by class lectures and non-computer oriented textbooks. 

Recently, this situation has changed due to the availability of time- 
sharing and telecommunication, facilities and to the cost rates being lower 
and accessibility much improved. More computer oriented textbooks in related 
fields like simulation [6,8], numerical analysis [1] and the systems approach 
[7] provide an initial base for the instructor in the development of computer 
aided instruction in the area of systems analysis. The aim of this paper is 



to descri'be such an attempt made retently in the Department of Agricultural 
Economics, Texas A&M University. 

The Course Objectives 
The objectives of the course were primarily to teach social science 
graduate students, mostly from economics and agricultural economics, the 
basic methodological and the quantitative tools'^of systems analysis and 
design [2,3,7j. Applied, as well as theoretical, consideration must be 
given so that planning and management of "real world" problems can be studied 
ar.^ solved. The student is expected to develop projects and homework assign- 
ments which involve simulation or to design a system to be checked and run on 
the computer along with a short essay type discussion summarizing the results, 
An example of such a homework assignment is provided in Appendix A. 

Course Content 

The content of the tv/o-semester graduate course is outlined as follows:' 
K Philosophy and methodology of systems science with references 

2. Mathematical modeling of dynamic deterministic systems with refer- 
ences [3,7,8]. ^ ' . 

3. Simulation techniques of dynamic-stochastic systems with references 
[7,9,10]. 

Control problems and optimization technique^ with references 

[1,5,7,8,9]. '/ 
The computer programs used FORTRAN IV ^ff^CSMP/360 [6] languages and were run 
on the Texas ASM University's 360/65 and/or 370/U5 computer^. 



The Computer Assistance in Force 

As mentioned above, emphasis was placed on having the students imple- 
ment segments oF systems science by using the computer. However, there was 
a wide range of backgrounds and experiences with regard to computer program- 
ming. This problem will probably persist at least for the next few years. 
Limiting the course attendance -o only those who have-had certain program- • 
ming experience, such as competence in FORTRAN, would have jeopardized the 
course due to lack of potential candidates quaHfied to take the course. 

In addition, a certain capability in computer programming and applica- 
tion is needed early in the game to meet the course objectives. To cope 
with this problem, a strategy was selected which will be discussed in three 
components: a) the simulation language, b) the FORTRAN extension and c) 
the system-algorithmic approach, 
a) The Simulation Language 

Careful consideration was given to the selection of the computer 
simulation language. It needs to have the capacity to handle dynamic and 
stochastic situations, and yet be easy to learn. The CSMP/360 language is 
well fitted for this purpose because It took only a few hours to teach the 
non-programmers enough features and structures of the language to allow 
them to advance on their own by assigning simple problems and access to the' 
User's Manual [6]. The capacity of this language to handle dynami c-stochas- 
tic systems is well established and documented [6, pp. 72-76]. The students 
were introduced to CSMP by having th*em generate wave signals using combina- 
tions of step, ramp, pulse and impulse functions. The students demonstrated 
their understanding of singular function concepts, as well as computer 
programming, in these exercises. The second step involved the solving of 
systems of n order differential equations by the various Laplace equivalent 



function blocks* This step also allov/ed non-linear interpolation to be 
* ' 2 

easily accomplished. The third step introduced stochastic proc^.sses which 
could be accomplished by the use of some simple FORTRAN programs. 

b) The FORTRAN Extension 

Due to the need to estimate model parameters, make data transformations, 
do matrix operations and generate stochastic variates distributed according 
to specific probability functions, some FORTRAN programming was required. 
A few FORTRAN help-sessions provided the students with the necessary level 
of skill. At this point, 'the students were taught to generate stochastic 
variates, Markov Chains and Monte-Car]o Simulations [9,10]. In each case, 
the students were furnished with the appropriate main computer programs 
which sometimes required the addition of certain subroutines in order to 
accomplish the assignment. The FORTRAN segments were also used for the 
simulation of discrete time models. The project assignments were carried 
out by small teams of students with at least one experienced in FORTRAN 
programming. The inexperienced were expected to learn enough FORTRAN so 
that the assignment could be carried out togethei". 

c) The System-Algorithmic Approach 

The complex project assignments were time consuming.to perform. There- 
fore, in order to reduce the time required by students and the computer cost 
of instrgction, an algorithmic approach was adapted whereby the overall 
program for each problem assignment was divided into program segments ab can 
be seen in Figure I . ' 
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Figure 1. 
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In this schematic, the solution for the problem requires three different 
algorithms to operate in a particular sequence dynamically. The goal is to 
train the students by using only the design of two specific algorithms* The 
overall program can then be segmented into a main program cotmecting three 
subroutines (macros or procedures in the case of CSMP) where the main program 
and the remaining subroutines are designed and applied by the students, who 
in turn devote a majority of their effort and time to the main assignment. 
Certain plot routines, estimations, optimizations or even complicated matrix 
operations may be beyond the capacity of the inexperienced programmer and 
consume too much of his limited time. This approach requires the instructor 
to carefully prepare these computer routines in advance. 

Some Observations 

This paper described briefly an effort to design a university level 
course in systems analysis with close contact and massive use of computer 
time in training the student along the entire spectrum of activities demand- 
ed by this discipline. In this case, the comput^i;^e comes the ^'laboratory 
room" for the systems scientist in much the same way students of chemistry or 
physics need their laboratory experience. On the average, the students 
spent a ratio of 2.5:1.0 in computer design programming, debugging and 
analysis versus class lecture time. It should be noted, though, that most of 
them did not take any otjher course and only shared their time with their on- 
going research activities. 

The results were encouraging and point to the need for further effort 
In the direction of additional computer exercises and the renK)val of much of 
the programming burden from the students' assignments by providing them with 
programmed segments. 
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Another significant step forward may be expected as computer satellite 
terminals become more available to students. Although both FORTRAN and 
CSMP are available through such terminals, the inclusion of interactive 
languages, such as APL, may provide additional convenience and flexibility 
in performing the course exercises and projects. 



Footnotes 



'^Assistant Professor,. Department of Agricultural Economics, Texas A&M 
Urii versi ty, 

Whe interested reader couId"get a detailed view of the course content 
by following the cited references or by contacting the author, 

2 . • 

See Appendix A for an example of a homework assignment which was 
typical at this level , 
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APPENDIX A 

t 

This appendix presents an example of a project assignment, requiring 
students to simulate a simplified case of pest population dynamics under 
control. The problem statement is given first, followed by a student's 
solution in terms of a block diagram, CSMP program list and pprtion of its 
output. 

Pest Control Simulation 

The Problem Statement 

Consider a single crop, single pest population (P(t)) interaction. The 
population dynamics is influenced by the fbl lowing characteristics: 

1. Initial infestation level P . 

o 

2. Natural mortal ity: Mj ' P(t) 

where Mo* = Pretreatment natural mortality rate 

Ml = Post 1st treatment natural mortality rate 

3. Eggs lay: E, * P(t) 

the emergence df adult pests occur with 20 days delay. 
, Adul ts migration 

the number of pests migrated into the field is a random variable 
normal ly distributed 

Ip = n(iJ,a") 

5. Number of pests killed due to treatment 

Kp = K^(x) • K^(k) • P(t) 

where K^Cx) = 1 - e "^^'^ 

K^(k) = f(k) (arbitrary function given below) 

K = 1,2,3 are days since treatment 
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6. Damage value Is proportional to P(t) and given by 

DMG - P(t) • D • B (per day) 
where D = per insect per day physical damage 
3 = product price 

7. Cost of treatment is given by 

COST t if X > 0 

{0 if X = 0 

where 

C = setup cost per treatment 

a = pesticide purchasing price plus direct application costs. 

Total accumulated cost = N(C+ax) where N is the number of treatments, 
ine goal is to develop a simulation system capable of making 
sensitivity analysis. 

Execute this Sssignmen't in two separate parts: 

Pa^^. Feasibility study including all steps up to programming (detailed 
block diagram included) . 

Part_B. After clearing your block diagram (a must), program it and run it 
with the following initial conditions and parameters: 

= 2500.0; P^ (population threshold for treatment) = 5000.0; = 0.125; 

= 0.100; E, = 0.15; M = 150.0; a = 75. 0; a = \.Ol\l\; X = 1.025; C = 3.0; 

3 = 0.7; a = 1.5; x = 3.8; f(l) = 0.7; f(2) = 0.5; f(3) = 0.3. 

let the final time = 100; DELT = 1.0; OUTDEL = KO 

Print at least the variables p(t) and total cost. Print-plot the variables^ 
P(t) and Kp. I^ethod of integration RECT. Maximize (approximately) the^profit 
by parameterization of the P^ and x variables. Make more of your own. 
assumptions if needed. 
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♦ ♦«*CPMTI NUnUS SYSTFM MODELING PI^OG«AM**** 

*■ 

VFRS ION' I • 3 **♦ 
TITLE PCST CONTROL SIMULATION 

i,^ >^«^X«* ******* 

INI TIAL 

*<c** ******* *x^* ************* ******* ********************** 

PARAMrTF.R PC = 250C • 0* NVC~ •125*M1=,1 ♦PC = 50 0C.0tX = (l .t4.»8. ) .f:l=.15 
CONSTANT B= . 7 ♦D=. DO I 1 C = 3.G » S-l'.f>f A = l c044 , W=l . 025 ^ 
* 

**<c *,{c«-J^**«*****«******** *****i!c** ********************* ***<'** *'X ******** ***^^ 

DYNAMIC 

*****>-*** **i?. *** V** *************** ***-**<^'i' ******** ***** A **************** ********* 
* 

PT=: I NTGRL (Pf I PP) 
RR=EL"KD-YY-KPCO-KP':)+ I P 
YY=PTj^Y 

Y=F CNSVi ( L , MO . MO i M I ) 
Yl =-CrJ''PAR( PT*PC) 
KP=Y1 '"KD*. /''^"PT 
Y2=X**W 

KD=1 .r-<rKP{-A*Y2) \ 

KPD1=Y1 *KD*.5*PT 

KPD = DC:LAY{ 1 , 1 .0* KPDl ) 

KPD ?- = Yl *. 3*Kn*PT 

KPDD=DtLAY( 2i2.C »<PD2) 

F2= ( Fl ) *(PT) 

EL^'^CL AY(20i:>0.0i-72) ^ 

LL- ( C+S*X ) *Vl ^2.0 

L=I NTGPL(0 • r ,LL) 

Z7 = Pr*DV3 

Z=INTGPL{0 .r. , ZZ) 

AC=74.L 

IP=GAUSS( 1 1 1 50. 1 75. ) 

**4t*^«* ************ ********* **)♦*********************************♦♦ ***********'^*^ 
TFPMINAL 

************** A* ******x^* ************************************************ ****** 

TIM^iP F INT IM=1 • D ♦ DLL T= 1 • i OUTO^L = ? • C 

* 

PPINT E? % KOtKP 

PPTHLT OTiAC f FL ♦ LL T Z7. % KP02 % KPOD ♦ P^ ♦ L ♦ IP ♦ KP 

MrTHOD TPAHZ 

fNO 

STDO 
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PROBLEM DURATION 0.0 



TO l.OCCOF 02 



VAPIA8LE 
PT 
AC 
EL 
LL 
22 

KP02 

KPOD 

RR 

L 

IP 



MINIMUM 
\.272\E 03 

0 • e 

O.D 
0.0 

8.9n50E:-01 

0.0 

O.C 
-1.8796F 03 

O.C 
-9. 24995 01 

0.0 



T IME 
1 . 9000E 0 I 
0 .0 
0.0 

0 .0 

1 .9000E 0 1 
0.0 

0.0' 

9.40GOE 01 

O.C 

0.0 

0 .0 



MAX IMUM 

5.0629E 03 

2.56C3E 02 

7.5493E 02 

9.000CE GO 

3.544CE OC 

9.8415E C2 

9.8416E 02 

5.4541F C? 

2.7000«=; 01 

2.304CF 02^ 

2.2964E 03 



TIME 
9.4000^ 01 
l.OOOOE 02 
9.9000E 01 
7.9000E 01 
9.40C0E 01 
9.4000F 01 
9.600DE 01 
9.9CC0E 01 
9.5000E 01 
7.9000E 01 
9.4000E 01 
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CO 
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MINIMUM PT VERSUS TIME MAXIMUM 

1.C865F 03 X = 8.D000F 00 5#0629E 03 

PT ' I -J 

2.5OC0E 03 + 

2.1694E 03 4- 
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